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Abstract 

Theoretical description of diffusion growth of a gas bubble after its nucleation in supersaturated liquid solution is 
presented. We study the influence of Laplace pressure on the bubble growth. We consider two different solubility 
laws: Henry's law, which is fulfilled for the systems where no gas molecules dissociation takes place and Sievert's 
law, which is fulfilled for the systems where gas molecules completely dissociate in the solvent into two parts. We 
show that the difference between Henry's and Sievert's laws for chemical equilibrium conditions causes the difference 
in bubble growth dynamics. Assuming that diffusion flux of dissolved gas molecules to the bubble is steady we 
obtain differential equations on bubble radius for both solubility laws. For the case of homogeneous nucleation of a 
bubble, which takes place at a significant pressure drop bubble dynamics equations for Henry's and Sievert's laws are 
solved analytically. For both solubility laws three characteristic stages of bubble growth are marked out. Intervals of 
bubble size change and time intervals of these stages are found. We also obtain conditions of diffusion flux steadiness 
corresponding to consecutive stages. The fulfillment of these conditions is discussed for the case of nucleation of water 
vapor bubbles in magmatic melts. For Sievert's law the analytical treatment of the problem of bubble dissolution in a 
pure solvent is also presented. 



1. Introduction 

This paper presents a theoretical description of diffusion growth of a gas bubble in liquid solution as a result of 
a considerable pressure drop (in the order of 10^ times). These conditions of bubble growth process are observed 
in magmatic melts during volcanic eruptions IT] |2l- After such a significant pressure drop the solution becomes 
strongly supersaturated; and homogeneous (fluctuational) nucleation of gas bubbles becomes possible. It is the growth 
dynamics of such bubbles that is the subject of the present paper. It has to be noted that the growth regularities of a 
solitary bubble are crucial for the description of the whole kinetics of phase transition in supersaturated solution ||3]|4l 

a. 

While describing gas bubble growth in supersaturated solution two rough approximations are traditionally made 

EiisiEiiiia: 

1 . The flux of the dissolved gas towards the bubble is assumed to be steady. 

2. The consideration is limited to bubbles with the radius that is large enough to neglect the Laplace pressure in 
comparison with the external pressure of the solution. 

In the present paper, following our preceding works ||8] 0, while exploiting the steady approximation, we take into 
account the Laplace pressure in the bubble considering the bubble from the very moment of its nucleation. Consider- 
ation of the time-dependent Laplace pressure in the bubble makes both gas density in the bubble and the equilibrium 
concentration of the dissolved gas at the surface of the bubble time-dependent as well. 
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Gas bubble growth in a solution taking into account the Laplace pressure was considered as early as in 1950 in 
the classical paper by Epstein and Plesset ifTOl . where the authors obtained the equation for the bubble radius as a 
function of time. In order to relate the equilibrium concentration of the dissolved gas and the solution pressure, IfTOl 
presupposed the fulfillment of Henry's law: i. e. the proportional dependence between these two values. Indeed, 
Henry's law is fulfilled for the solution of CO2 in magmatic melt Q; however, it is not valid for the solution of 
H2O vapor 13] |6l. In this case, which is crucial for practical reasons, Sievert's law is observed: the equilibrium 
concentration of the dissolved gas is proportional to the square root of the solution pressure IITTI . Sometimes both 
cases are referred to as "Henry's law", but in the present paper, in order to avoid confusion, we will use the term 
"Sievert's law" for the case with a square root, following e. g. lfT2l . Here we analyze the bubble growth dynamics in 
both cases: for Henry's and Sievert's laws. For Sievert's law we obtain the equation for the bubble radius as a function 
of time analogous to IfTOl . which has not been obtained previously. 

Our analysis shows that, irrespective of the law applied to gas solubility, three characteristic stages can be marked 
out in the growth dynamics. During the first stage the bubble radius is growing with an increasing rate. On the second 
stage the growth rate decreases. The third stage, when the growth rate continues to decrease, begins when the Laplace 
pressure inside the bubble becomes comparable with the external pressure of the solution. We demonstrate that during 
the first two stages the time dependence of the bubble radius is different for the cases when either Henry's and Sievert's 
laws are fulfilled, while during the third stage this distinction is no longer observed. 

For both Henry's and Sievert's laws we obtain intervals within which the bubble radius changes on each stage, as 
well as time limits and conditions when the steady approximation is applicable. We show that, as the radius of the 
bubble increases, the steady condition becomes stricter; and, consequently, as a rule, the steady regime of a multistage 
bubble growth gradually gives way to the nonsteady one. We obtain analytical expression for the time when Laplace 
pressure influences on bubble growth vanishes and therefore substantiate the estimation of this time made in |[T3l . 
After this time passed the bubble growth reaches a self-similar regime lfT4l . We also analyze whether the steady 
approximation is applicable to the case of gas bubbles in magma described by Navon ||6l and Chernov et al. ||3] for 
large radius of a bubble (neglecting Laplace pressure). Besides that, we present the analytical description of bubble 
dissolution in the pure solvent for the Sievert's solubility law, which was not presented before in literature. 



2. Equilibrium Concentration of the Dissolved Gas 

Let us consider a gas dissolved in a liquid. The solution was initially in an equilibrium state at temperature T 
and pressure Pq. The concentration of the dissolved gas in the solution under such conditions will be denoted as no- 
Then we instantly relieve the external pressure to value 11 in such way that solution becomes supersaturated. The 
temperature and volume of the solution remain the same, thus value no still serves for the dissolved gas concentration. 

It is more convenient to express the state of the solution in terms of dimensionless variables: supersaturation f and 
gas solubility s defined here via 

(1) 

noo 

. = -f^, (2) 

where nco is the equilibrium concentration of dissolved gas at the external pressure 11, k is the Boltzmann's constant. 

When some time passes after the pressure drop, a gas bubble nucleates and begins growing regularly. Following 
|[3) we assume that the bubble is in mechanical equilibrium with the solution, and its dynamics is governed only by 
diffusion process. This assumption will be discussed in Appendix A. The radius of the bubble will be denoted as R. 
We consider the situation when the solvent is in its stable liquid state; therefore, we assume that the bubble consists 
of gas only, but not of the solvent vapor The gas in the bubble is considered to be ideal. 

When the bubble is studied after some time fo since its nucleation, its radius complies with the strong inequality 

R » 2o-/n , (3) 

where cr is the surface tension of the pure solvent (it is true while the solution is considered as diluted). Eq. Q allows 
us to neglect the influence of Laplace pressure on the bubble growth. Therefore, we can write the following equations 
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for the pressure in the bubble Pr and for the equiHbrium solution concentration near the surface of the bubble riR. 



P« = n , (4) 

riR-rioo ■ (5) 

The subscripts oo denotes that the equilibrium concentration iioo is related with the equilibrium near the flat surface of 
phase separation {R — > oo). From Eq. (HJi it follows that the gas concentration in the bubble Ug is constant. Using the 
ideal gas law we have 

n„ - — . (6) 

^ kT ^ ^ 

When Eqs. (|4]i, and (|5]l are fulfilled and, therefore, Laplace pressure is negligible, the bubble dynamics is evident 
for the case of steady-state diff'usion and can be even described analytically for the non-steady case |[T3l[T4l . 

From the moment of bubble nucleation and till Eq. Q becomes valid, Laplace pressure influences bubble growth. 
Thus both quantities Pr and «r become radius-dependent (and, therefore, time-dependent). For Pr now we have 

PR-n + ^. (7) 

In order to write the equation for we need to know the solubility law. For the simplest case of Henry's law the 
equilibrium concentration of dissolved gas is proportional to the corresponding pressure 

^ = (8) 

However, Henry's law is fulfilled only in such systems where there is no molecules dissociation during gas dissolution. 
In another important case when the gas molecules completely dissociate in the solvent into two parts the Henry's law 
is replaced with so-called Sievert's law (see e. g. ifTSi ). and, therefore, Eq. ^ is replaced with the following one: 

(9) 

Sievert's law is fulfilled for water vapor dissolved in a silicate melt IfTTI . Such solutions are important both for glass 
production lfT2l and for volcanic systems ||6]|3|. 

The replacement of Eq. (O with Eq. (|9]l means that the boundary condition for the gas diffusion problem will be 
different for Henry's and Sievert's laws. The change of boundary condition, as we will see further, leads to the change 
of bubble dynamics. 




3. Bubble Dynamics Equation 

After the nucleation of the bubble, when its growth can be considered regular (i.e. the bubble cannot be dissolved 
by fluctuations), its growth is governed by the diffusion ffux of gas molecules into it. In this paper we will study 
the case when the diffusion flux can be assumed as steady. The conditions when such approximation is valid will be 
discussed further. 

Taking into account the equality of gas concentration at the bubble surface to the equilibrium concentration hr and 
the equality of gas concentration far from the bubble to the initial concentration no, we can write a simple expression 
for the steady diffusion flux density jo 

no-riR 

JD = D — - — . (10) 

K 

Here D is the difliision coefficient of gas molecules in the pure solvent (we assume that the solution is diluted). 

Now let us write the expression for the number of gas molecules in the bubble. Exploiting the ideal gas law and 
using Eq. dTji we have 

20-1 



47r , 
3kT 



n + 

R 



(11) 



3 



Differentiating Eq. (fTTb . using Eq. (|6|l, we obtain 



dN jdR 



dt 



dt 



1.^ 
R 



where 



is the characteristic size of the bubble. 

Material balance between the dissolved gas and the gas in the growing bubble gives us the following equality 

dN , 

Substituting expression for diffusion flux density Eq. ( fTOl i and the rate of change of the number of molecules Eq. (fT2l i 
into material balance equation Eq. (fT4l i we have 



(12) 
(13) 

(14) 



dR 
'~dt 



1.^ 
R 



= D 



n o - riR 
R 



(15) 



Or, exploiting Eqs. (|2]l and (|6]l, we equivalently have 

RR 



1.^ 
R 



= Ds 



no - Wfi 



Ho 



(16) 



Using Eqs. ([TJ, d?) and, correspondingly, dS) and the fraction in the r h. s. of Eq. ( fTSI l can be expressed as 



ng-riR _ 2cr 

«oc ~^ 



for Henry's law, and 



no - n« ^ , /, 2o- 
= f + 1 - A/ 1 + 7^ 

Moo V n/? 



for Sievert's law. After exploiting in Eqs. (ITTt or (IT¥t the definition of (Eq. (fT3ll). we have, correspondingly, 

no - Mr 



3^ 
^ 2 R 



and 



no - nfi 



^+1 



'iRo 



noo y 2 R 

Eqs. (fT9] ) and ( l20b allow us to rewrite the equation of bubble dynamics Eq. ( fTSI l in the final form, namely 



RR 



1.^ 
R 



Ds 



3^ 

^ 2 



for Henry's solubility law; and 



RR 



1.^ 
R 



<r + 1 - a/i + 

^ ^' 2 R 



(17) 



(18) 



(19) 



(20) 



(21) 



(22) 



for Sievert's solubility law. 

While Eq. ^T\\ was obtained as early as in 1950 in paper by Epstein and Plesset IfTOl . Eq. ( l22b was not obtained 
previously. In paper by Cable and Frade it was presented only a special case of Eq. ( l22l i. when ( - - \: 



RR 



1.^ 
R 



2 R 



(23) 



(Eq. (34) in lfT2ll ). Such value of supersaturation comsponds to the dissolution of gas bubble in the pure solvent. 
And even for this special case of Eq. (l22l i analytical solution was not obtained. We present the analytical solution of 
Eq. (|23] | in Appendix B. 



4. Critical Bubble and Initial Conditions for the Bubble Growth 



The obtained equations for the bubble growth dynamics Eqs. (|2TI ) and ( l22b can be applied to the cases of both 
homogeneous and heterogeneous nucleation. Below we will consider only the homogeneous nucleation case. 

Homogeneous nucleation of a gas bubble in a supersaturated solution means fluctuational mechanism of its ap- 
pearance and needs a significant pressure drop (in the order of 10^ times). It is these conditions of bubble nucleation 
that take place in magmatic melts during volcanic eruptions 13) • Since we decided to describe only regular growth of 
a bubble, we need to exclude the very process of nucleation from our examination and consider a bubble only when it 
is already supercritical. 

Under the notion of critical bubble we understand, as usually ||4l[T6l, such a bubble which radius Rc corresponds 
to the extremum of work of bubble formation. Critical bubble is in mechanical equilibrium with solution at the initial 
pressure Pq and in chemical equilibrium with solution with concentration no 

Pr, = ^0 , (24) 

wr, = no ■ (25) 

These two conditions together with the solubility law unambiguously define the value of Rc- Condition ( l24b . using 
Eq. Q, evidently gives us 

la- 

Rc = , (26) 

Po-n 

and condition ( IZST i allows us to connect value Pq with the supersaturation 

Substituting Eqs. (l24l i and (IZST i in Eq. (jSj (Henry's law) and Eq. (|9j (Sievert's law) we have correspondingly 



»o _ Pq 

rioo n 

and 

no _ [Pq 

Hoc \ n 

For Henry's law we have evidently from Eq. (l27t using Eq. ([1} 



(27) 



(28) 



Po = n(^+l). (29) 

And, substituting Eq. ( |29] l in Eq. ( |26] |. we have 

2cr 

Rc^ — , (30) 

the well-known expression for the radius of a critical bubble (e. g. Bl fTSl ). 
For Sievert's law, when Eq. (l27t is replaced by Eq. (l28t . we have 

Po-nM\ (31) 



or, using Eq. ([TJ, 



Ho 



Po = n(^+l)2. (32) 



Finally, substituting Eq. ( |32] | in Eq. {TM for the critical bubble radius, we have 

2cr 

~ n[(^ + 1)2-1] • ^^^^ 

Using Eq. ( flj] ) and Eq. ( |30] | or Eq. (|33] | correspondingly we can write the relations between R^- and Rc for both 
Henry's and Sievert's laws 

Ra = M , (34) 
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R^^^[(^+lf-l]Rc, (35) 

which will be exploited later. 

Here it is important to make a remark regarding the quantity s, the gas solubility. Eqs. ([TJ and (|2]i give us 

npkT 1 

■s - 7 • (36) 

n ^+ 1 



If Henry's law is fulfilled, we have ( + I - Pq/H and, therefore, 

nokT 
s - , 

^0 

s is a tabular value defined only by the initial state of the solution (before the pressure drop). 

Since we consider Sievert's law, we have (see Eq. (l32t ) ^ + 1 = yjPolYl; thus instead of Eq. ( l37T i. we obtain 

nokT 



(37) 



(38) 



Here it is convenient to introduce constant K as the coefficient of proportionality in Sievert's law: K = no/ VA)- It 
allows us to rewrite Eq. (l3Ft in the following form: 

kT 

s = K— . (39) 

Eq. ( [39] l shows us that in the case of Sievert's law the solubility value s depends on the final state of the solution. In 
the current paper a solubility value for Sievert's law is understood as gas solubility at the final pressure 11 i.e. after the 
pressure drop. 

To deal with the dynamics equations (ISTT l and (l22l l one needs to provide it with a reasonable initial condition, i. e. 
to choose the value of constant initial radius of the bubble, in the following equality 

Rit)\,=o = Ri ■ (40) 

In papers IfTOl [T2l there was no special meaning assigned to the value of /?,: the reason of appearance of the bubble 
was totally excluded from discussion. Here we assume that the bubble nucleated fluctuationally, i. e. crossed the 
barrier corresponding to radius R^. It means that has to be not less than R^, but evidently we cannot use the radius 
of a critical bubble R^ as the initial value for the radius. 

A bubble that nucleates fluctuationally in the solution is capable of regular growth if it passed and, moreover, 
moved away from the near-critical region where fluctuations are still strong enough. Thus we choose, following JS], 
Ri = 2Rc, i. e. 

i?(f)l,=o = 2R, . (41) 

Such value guarantees the absence of fluctuations and, as we will see further, provides us with convenient expressions. 

The necessary condition for the fluctuational nucleation of a bubble is a high pressure drop (Po/Il ~ 10^) and, 
consequently, high supersaturation ^ ~ lO'' for Henry's law and ^ ~ 30 40 for Sievert's law. Further we will use the 
following strong inequalities for the supersaturation: 

^ » 10 (42) 

for Henry's law, and 

^ » 1 (43) 

for Sievert's law. 
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5. Three Stages of Bubble Growth 



Before we find the explicit solution of equations (ITTT i and (l22l i. let us qualitatively describe the change of the 
character of bubble growth process with the increase of its size. This will allow us to mark out the representative 
stages of growth and to determine corresponding characteristic bubble sizes. Duration of the consecutive stages and 
the character of bubble radius time dependence on each stage will be considered in the next section. It has to be 
noticed that the stages of our interest do not have anything in common with the stages of evolution of the whole 
ensemble of bubbles during the decomposition of liquid solution supersaturated with gas. 

It will be more convenient to consider both cases of solubility laws separately. 

5.1. Henry's Law 

Let us rewrite Eq. (ISTT i in the equivalent form which will be more appropriate for its analysis: 

* = °<('-|)^(tt^)- 

where we took into account Eqs. (flJT l and (l30l l. 

Each of the three co-factors dependent on R emphasized in the right hand side of Eq. (l44l i describes its physically 
different contribution to the dynamics of the supercritical bubble growth process. Co-factor 1 - Rc/R, increasing 
with R, corresponds to the fast (the scale of change of R is R^) increase of the driving force of the process (the value 
nQ-nft) with the growth of R. Co-factor IjR, decreasing with the growth of R, describes, as it is seen from Eq. (fTOt . the 
contribution related to the decrease of the gradient of solution concentration near the bubble surface, which decreases 
the bubble growth rate with the growth of R. Finally, co-factor 1/(1 + R,t/R), which increases with the growth of R, 
takes into account the counteraction of Laplace pressure to the bubble growth, i. e. the fact that the bubble growth is 
facilitated by the reduction of Laplace pressure with the growth of R, while other factors are equal (the scale of its 
change is Ra-)- Notwithstanding the mentioned reduction of counteraction, the resulting contribution of the last two 
factors always leads to the deceleration of growth with the increase of bubble size. 

From Eqs. ( l34l i and (l42l i it can be seen that R^r » Rc- Using this inequality, Eq. (l44l i and boundary condition (HTb . 
we have that the growth rate of the bubble radius in its dependence on this radius has to reach the maximum value 
achieved at certain bubble radius from the interval 2Rc < R < Ra-- Thus it is natural to consider the growth in the 
following interval of sizes 

2Rc<R< R^ (45) 

as a first stage of bubble evolution, where the determining factor is the increase of the driving force of the growth. At 
this stage bubble growth goes with the increasing in time rate, reaching its maximum at R„,. 

In order to obtain R„i we will consider the rate of bubble growth as a function of its radius and differentiate both 
parts of Eq. (l44l i with respect to R: 

dk Rc (Rc + Ra) -(R- Rcf 

— = DsC-^^ — — — ■ (46) 

dR ^ RHR + Ra? 



The quantity Rm is defined by the extremal condition 

dk 

dR 

which using Eq. (|46] | leads to the following result 



= , (47) 

R=R„ 



R„^Rc + [Rl+RcRa)"\ (48) 
Taking into account the strong inequality Rcr » Rc, we can simplify the obtained expression for/?,,,: 

Rm - (RcRa)''^ . (49) 

The second stage of the process will be when the bubble growth occurs within the interval of sizes 

R„,<R< R,r . (50) 
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During all this stage, as R„ ^ (/?c/?o-)''^ » Re, it is already vaUd that 7? » Rc, thus the driving force no - «/? remains 
practically constant. As a result, bubble growth decelerates, although, as it was noted above, the counteraction of 
Laplace pressure to the growth gradually is attenuated. The Laplace contribution Icr/R to the pressure inside the 
bubble decreases during the second stage by {Ra-jRcY^'^ times and at the completion of this stage becomes comparable 
with the external pressure EL contribution (from Eq. (fTST i we have strict equality la-jRcr - 311/2). It is this physical 
condition that defines the completion of the second stage. 

On the subsequent, third stage, which corresponds to the interval of sizes 

R>R.r, (51) 



monotonous decelerated bubble growth continues. At the same time, the role of Laplace pressure continues to decrease 
gradually, and the pressure inside the bubble approaches to a constant value equal to the external pressure 11. As it 
will be shown further, this process is rather protracted, so the concluding phase of the third stage, when the pressure 
inside the bubble practically does not change and the application of self-similar solution |[T3lfT4l is possible, comes 
only in the interval of sufficiently large sizes of the bubble, when the condition R » Rg- is satisfied with a certain 
reserve. 



5.2. Sievert's Law 



Let us investigate behavior of bubble growth rate R with the increase of bubble radius for Sievert's solubility law. 
From Eq. (l22l l evidently stems. 

Ds 

^+1 



R 



R + Ra 



3Ro 



2 R 



(52) 



Eq. ( l52b can be also written in the form similar to Eq. (04), to mark out three co-factors. The character of bubble 
growth rate is presented graphically in Fig. [1] 
Diff'erentiating Eq. ( |52] | by R we obtain 



Ds 



dR 



dR (R + Ra-Y 



3 RAR+Rg) 
4 



R2 



^+1- 



2 R 



(53) 



It can be easily seen, that, as it is for Henry's law, here bubble radius growth rate as a function of variable R also has 
the only maximum. Denote the corresponding radius as R^, we can obtain its value from Eq. (l47l l. Using Eq. (l53l l we 
rewrite this equation as 

4 



RR,r +Ri + -R'll + -—]- -R 
2 R 



1 + 2^(^+1)^0 
2 R ^ 



Assuming that the sought quantity R„ is considerably less than Rg-, we will use for R in Eq. 
inequality 3Ra-/2R » I. Below we will need this inequality to be even more stronger 



(54) 
the strong 



R^/3R » 1 . (55) 

It will allow us to omit the second term in the brackets and the second addend in the square root in Eq. (l54l . So this 
equation can be rewritten as: 

^ V^(^ + 1 )R^'^ - 3RR^ -Rl^O. (56) 

Using Eq. ( fSST i we can also omit the second addend in Eq. ( |56] | in comparison with the third one. After such a 
simplification Eq. ( |56] | becomes solvable; and for R„, we have: 



Rm — 



</3 R^ 

2 ({+ 1)2/3 • 



(57) 
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Substituting Eq. dSTb in the second addend in Eq. d56l t we can find a first-order correction to R,^, in Eq. dSTl ). We have 



Rm — 



V3 Ra 



2 (^+ 1)2/3 



1 + 



3<f3 



2 (^+ 1)2/3 j 



2/3 



(58) 



Numerical solution of Eq. ( |54] | for ^ = 30 gives for inaccuracy of approximate solution Eq. (ISFt the value less than 
2%, and this value evidently decreases with the increase of 

As soon as we explained the behavior of the value R and found the value of R,,, we can use the ideas proposed for 
Henry's law to determine the stages of bubble growth. On each stage the character of bubble growth is different from 
any other. The three consecutive stages of growth are still defined by Eqs. i45[ . dSOl l and ( BTT l. 

In the next two sections we will obtain the time dependence of bubble radius R for each stage and duration of each 
stage. Also we will obtain conditions of steadiness of bubble growth for each stage. 



6. Time Dependence of Bubble Radius 

Let us now solve the differential equations (ISTT l and (l22l i for the time dependence of the bubble radius with the 
initial condition Eq. ( l40t of homogeneous nucleation of the bubble. As it was in the previous section it is more 
convenient here to consider Henry's and Sievert's laws separately. 

6.1. Henry's Law 

At first we will rewrite Eq. (l2Tt in the form which is appropriate for integration: 

RR + {R^ + R,)R + (R^ + Rc)Rc-^^i^ = Ds( . (59) 

R/Rc - i 

Integrating Eq. f59[ . we obtain 

y +(/?,, + R,)R + (R^ + R,)Rc In - 1 j = D< (f + t) , (60) 

where r is the constant which has time dimensionality and which is defined by the initial value of radius at time f = 0. 
Using initial condition (|40|) . from Eq. (l60l l we find 

2RARa + 2Rc) 

-57^ ■ ^'^^ 



Excluding time r from Eq. (l60] t by means of Eq. ( 1611 , we obtain 



47^^ / 

" +{Rcr+ Rc){Rcr - 2Rc) + {Ra + Rc)Rc In 11= DsC,t . (62) 



\K(- 

Validity of the general relation Eq. (l62T l. which strictly takes into account Laplace pressure influence on the bubble 
growth process, is limited only by the condition of applicability of steady approximation Eq. ( fTOl l for the diffusion 
flux. Equation ( |62] | does not imply the smallness of quantity Rc/Rcr which follows from Eqs. (|34] | and (l42l i. Eq. ( |62] | 
complies with the results obtained in papers ifTOl [TSI for the particular case of steady growth of a bubble. 

Now let us consider the third stage of bubble growth, when R > R^. First of all, let us note that at the end of the 
first stage, when the bubble radius R approaches the value (RcRa-y^^, and with even more assurance on the second and 
the third stages, one can neglect the logarithmic addend in Eq. ( l62b . Moreover, since during the third stage the main 
contribution in the 1. h. s. of Eq. (|62] | gradually tends to the first addend, equal to R^/2; and the contribution of the 
second addend (influence of Laplace pressure) decreases, neglect of logarithmic contribution becomes fairly justified. 
As a result, Eq. (l62l l conformably to the third stage of bubble growth can be written in the form of 

R^ 

— + RcrR = Ds^t . (63) 
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6.2. Sieve rt's Law 

Unlike the case of Henry's law (Eq. (ISTT i) solution of Eq. (l22T l in general case is too cumbersome and, as we will 
ensure in the current section, when inequality (l47t is fulfilled, is even not necessary. Here we present solutions of this 
equation for two particular cases; 

R « R^ (64) 

- the first and the second stages and 

R»R, (65) 

- the second and the third stages. When strong inequality Eq. (l43T l is fulfilled, these two cases cover the whole range 
R > 2Rc of regular growth of the bubble radius, and that is the reason why general solution of Eq. (l22l i is not necessary 
for the system under consideration. 

Let us begin with the case when inequality Eq. ( |64] | is fulfilled. In this case we can omit 1 in comparison with 
Rcr/R in the 1. h. s. of Eq. ( l22l i and we can also omit 1 in comparison with the fraction 3Ra-/2R in the i. h. s. of this 
equation. Thus we have 



. Ds 
R^ — 



; 3 



(66) 



Separating variables and exploiting Eq. (|35> for Rcr under the square root, we can rewrite Eq. (I66t in the form which 
allows its integration: 

= (67) 
1 - ^/RJR R<r 

Integrating Eq. ( |67] | with initial condition (l40l i. we finally have 



R-2Rc+Rcln^^ - lj + 2V^(V^- ^^IR^+RM 



f V2 + 1 V^- 

V2- 1 v^. 



Ds 

= — (^+l)f. (68) 



This formula is different from Eq. (|62] |. which means that, when inequality (|64] | is fulfilled, there is a significant 
difference in the character of growth between Sievert's law and Henry's law. 

Now let us proceed with the other case. At first, using Eq. d35t we can rewrite Eq. (|22j equivalently in the form of 

=D.(^+1). (69) 



1 _ /^_ + «£ 

This form makes it obvious that, when strong inequality (l65T l together with inequality (l43l l are fulfilled, the whole 
square root in the denominator of the 1. h. s. of Eq. (l69t can be omitted in comparison with 1. Therefore, we have 



R\R + RA^Ds{l^+Y) . (70) 

This expression can be easily integrated. But the use of the initial condition Eq. (|40] | is not just as a result of the 
fulfillment of strong inequality ( l65l l. There is arbitrariness in the choice of the initial condition for integration of 
Eq. (iTOl l. but the most convenient is to choose a condition at such an "average" radius which simultaneously satisfies 
Eqs. (l64l i and ( |65] ). e. g. 

/J(f)U,, ^Ra= ^^RcRa ■ (71) 

Due to Ra <*; Rcr, we can use Eq. (l68b to obtain the explicit value for time ta defined in Eq. dTTl i. Using inequality 
Ra » Re, Eq. (I68l l gives us 

RgR. „2) 

Now, integrating Eq. (ITOl t with initial condition dTTT t, we have 

-—^+RAR-Ra)^Ds{i+\){t-Q, (73) 
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Figure 2: Growth rate of bubble radius squared dP?- jdt (measured in dimensionless units 2Ds() as a function of R. Solid curve con'esponds to 
Henry's law, Eq. (44). Dashed curve corresponds to Sievert's law, Eq, )52t . For both solubility laws Po/n = 10^. 
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or, rewriting Ra using Eq. ( ITTI ) and tg using Eq. ( l72l i. we have 



^ (74) 

With the increase of R the contribution of the third addend in the 1. h. s. of Eq. (l74l l decreases and alR > R,y (on 
the third stage) we can akeady write 

^2 

— + R,^R ^ Ds{^ + l)t . (75) 

If we also used in Eq. ( fTST l strong inequahty ( l43T l, Eq. (iTST i will become identical to Eq. (l63T l. It means that, when 
inequality R > Rg- is fulfilled, any difference in the character of growth between Sievert's law and Henry's law 
disappears. 

When the bubble radius becomes as large as 

R:»R,r, (76) 
Eq. ( l63T l transforms into the well-known Scriven's |[T4l self-similar dependence for the steady-state case 

R^ = IDs^t , (77) 

when dR^/dt = const. This trend of Eq. ( l63T l toward Eq. (ITTI i was discussed in detail in ||8]. This trend is presented 
graphically in Fig. |2] 



7. Duration of the Consecutive Stages 

7.1. Henry's Law 

Eq. ( |62] |. giving the explicit dependence of bubble radius on time, allows us, in particular, to find characteristic 
times corresponding to consecutive stages. According to Eqs. ( l45T l and (l48T l. at the end of the first stage, the bubble 
radius reaches the value R,„ = {RcRo-Y^^. Substituting value R - R„, to Eq. (|62] | and considering that (by virtue of 
inequalities Rc <sc R,,, <s: R,j) the main contribution to the 1. h. s. of Eq. ( |62] | is made by the second addend, we obtain 
the expression for the first stage duration ti 



Ds(\Ra 

Using Eq. (l34l l. expression dTST l can be also presented in the form 

3\'^' Rl 



(78) 



As it follows from Eq. (|79]l, with the increase of initial supersaturation of the solution, the first stage duration decreases 
proportionally to 1 f^^^^. 

The second stage of bubble growth starts at the time point ti and finishes at the time point t2 defined by the 
condition = Ra-- Substituting value R - R^^ into Eq. ( |62] | and considering that, by virtue of inequality R^ Rg-, 
the main contribution to the 1. h. s. of Eq. (|62] | is made by the first and the second addends, we obtain the expression 
for t2 

3 Rl 

f2 = -— . (80) 

As one can see from this expression, f2 dependence on initial solution supersaturation is defined by multiplier 1 
The duration of the second stage is much longer than the duration of the first stage, as from Eqs. (|79] | and (IHOt it 
follows that 



t2 IH''^' 



h \2 



» 1 . (81) 



Now let us consider the third stage of bubble growth, when R > R,j and f > f2- First of all, let us note that at 
the end of the first stage, when the bubble radius R approaches the value {RcRa-Y^^, and with even more assurance on 
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the second and the third stages, one can neglect the logarithmic addend in Eq. ( |62] | (it was taken into account earlier, 
when Eqs. ( fTST l and ( l80l l were obtained). 

Due to Eq. ( BTI ) the duration of the third stage of bubble growth is infinite. But it is reasonable to estimate the time 
from the beginning of the third stage and to the moment when bubble radius reaches the value defined in ITSl . 

Ro = 20 X 2cr/n . (82) 

It is assumed a priori in |[T3l . that when radius reaches value /?(> the influence of Laplace pressure on the bubble 
growth is negligible. 

As it was found in the previous section, on the third stage of bubble growth the time dependence of bubble radius 
is given by the simple equation, Eq. ( l63T l. Let us introduce time tj, as a duration of bubble growth in the size interval 
Ra- < R < Ro- Using Eqs. (l63T l and (|82] | we evidently have 

f3 480— . (83) 

The duration of the third stage is much longer than the duration of the second stage (and, moreover, the first stage), as 
from Eqs. (jSOll and (l83]l it follows that 

- ^ 320 . (84) 

Eq. ( l84l l and strong inequality (ISTT i allow us to evaluate the whole duration of bubble growth in the interval of sizes 
2Rc < R < Ro as t^,. In |[T3ll the estimation of this duration fo was given by the following equation 

fo ^ — — ■ (85) 
2Ds( 



Using Eqs. ( fT3T l and ( l82l i we have to ^ 450;^; and therefore 

6% . (86) 



h - to 



7.2. Sieve rt's Law 

Using Eqs. ( l69b and ( fSST i we can obtain expressions for duration of the first two stages fi and f2 for the Sievert's 
law analogous to Eqs. ( iTSl l and dSOl l. These expression are cumbersome, while qualitatively the time-scale hierarchy 
for Sievert's law is similar to the one for Henry's law. Therefore we will not present these expressions here. 

As long as on the third stage the dynamic equation of bubble growth for Henry's law Eq. (l67t and for Sievert's 
law Eq. (l75b are exactly the same (we need to account f » 1 also), all the results presented above for the third stage 
duration are vahd for Sievert's law also. 



8. Steady Flux Conditions 

The diffusion flux of molecules toward a growing bubble can be considered steady when the bubble growth is slow 
enough in comparison with the "diffusion cloud" growth. To be more exact, the radius of the bubble has to be much 
smaller than the radius of this cloud, the diffusion length. We can express it as 

R « (DtRf^ , (87) 

where is the characteristic time of the bubble radius change, = R/R the time in which the bubble radius changes 
significantly. Evidently, Eq. dSTI i can be rewritten as 

[RRIof^ « 1 . (88) 
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We can make this condition more explicit by means of Eq. (ISTT i and Eq. ( l22b that give us correspondingly: 



for Henry's law, and 

f I — Tr\"^ 



1 + ^ 



«; 1 (90) 



for Sie vert's law. 

Now exploiting Eqs. ( |89] l and ( |90l l let us obtain the conditions for diffusion flux to be steady on each stage defined 
above. We will write these conditions as inequalities for the value of solubility, not for supersaturation. The value of 
supersaturation is already fixed by Eqs. (l42t and (l4Jt . 

Obviously the larger the bubble is the more strict condition for steadiness is (see Fig. |3]l. The general condition 
for bubble growth to be steady at any time is the condition at /? » R,r. From both Eqs. (|89] | and (|90] | we have 

s'" « (I)''' . (91) 

This condition is sufficient for steadiness on the third stage for both Henry's and Sievert's laws. 

For bubble growth to be steady during the whole second stage it is sufficient to be steady atR = R^r- This condition, 
using Eq. (l42l i. transforms Eq. (l89T l 

s"' « (^) . (92) 

Exactly the same result will be for Sievert's law (see Eq. (1431) and Eq. ( |90l l ?AR - Rg-). 

For bubble growth to be steady during the whole first stage it is sufficient to be steady at R = R,„. It is to be 
reminded that the quantity R,„ has different values for Henry's and Sievert's laws (see Eq. (29) and Eq. (fSTll). For 
Henry's law, using Eq. (l35l l this condition leads to 

^'"«iir- (93) 



3( 



And for Sievert's law to 



2 



^"'«[Yf- (94) 

Obtaining Eq. ( |94] | we used not only Eq. ( |43] ), but even more strict condition ^^'^ » 1 . This condition is still fulfilled 
than ^ ~ 30 H- 40. 

Finally, let us also write the condition of steadiness at the very beginning of the regular growth of the bubble, 
exploiting Eqs. ( [89] l and ( |90l l with R = 2Rc. For Henry's law we have 

s''^«ll\'\ (95) 



and for Sievert's law 

^1/2 ^ ^1/2 _ 

While deriving Eq. ( |96] l. the numerical coefficient ^3^1-1/ V2)) ^ in its 1. h. s. was replaced with 1 for shortness. 

Let us mention the following interesting observation. Since for R » Ra-, the bubble dynamics is exactly the same 
for both Henry's and Sievert's solubility laws, the condition of steady growth for R » R,ris also the same (see Eq. Wi} 
above). For the case of homogeneous nucleation, when the pressure drop Po/Il ~ 10^, the steady condition at /? » /?o- 
as a rule is violated in cases of both Henry's and Sievert's laws. For Henry's law, when s ~ 10"^, it is violated due to 
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high supersaturation values ^ ~ 10^. For Sievert's law, when corresponding supersaturation values are significantly 
less ^ ~ 30 H- 40, values of gas solubility are significantly higher than for Henry's law (see next section), and that is 
the reason of the violation of the steady condition. While for Henry's law solubility 5 is a tabular value (see Eq. (iJTTi). 
for Sievert's law it can be adjusted via settlement of the final pressure 11 value. It is evident that condition ( |9TI ) can 
be satisfied at the given solubility value when one decreases the solution supersaturation ^ (the case of heterogeneous 
nucleation). 

In order to satisfy the condition ( |9ll of steady growth the value of gas solubility s has to be decreased. From 
Eq. (l32l i for the supersaturation we have 

^=-^-1. (97) 
Than, using Eqs. (|39T l, (|97T i. let us rewrite product s( in the following form 




(98) 



Eq. ( |98] l shows us that in order to weaken limitation Eq. (|9Tt one needs to increase the final pressure H, leaving the 
ratio /"o/n constant. 



9. Bubble Growth in Volcanic Systems 

This section contains the analysis of the steady growth condition obtained in the previous section for bubble 
nucleation in volcanic systems. Previously, e. g. in papers Q and ||6l the study of such systems exploited steady 
approximation without the analysis on its applicability. 

For our evaluations we will use parameters from (|3] for the case of homogeneous nucleation of water vapor 
bubbles in a magmatic melt. We have 

Po = 100 MPa n = 0.1 MPa T = 1150 K w = 3% p„, = 2300 kg/ml (99) 

Here w is gas mass fraction of the dissolved gas (water vapor), and is magma density. 

Let us express the values of s and ^ using data given. From Eq. ( |97] l, in accordance with Eq. ( |99t , we have ^^31. 
Then we need to calculate no and substitute it into Eq. ( l38l l to obtain the value of solubility. As long as we are given 
mass density of magma and mass fraction of gas, it is convenient to write 

«o=Po— , (100) 

where po is mass density of the dissolved gas, A^^i =6x10^^ moP' is the Avogadro constant and fi - 1 .8 x 10"^ kg/mol 
is the molar mass of the dissolved gas (water). Finally, we need to express the mass density of the dissolved gas. 
Evidently, we have 

Po = wp„ , (101) 

and, therefore, 

NAkT 

s = wp„ — —= . (102) 

Using data given by Wq. ( |99] l in Eq. ( I102l i. we have s ^ 12. 

Now we can see that for volcanic systems, where the pressure drop is of the order of 10^ and solubility is more 
than 1, both conditions (|9TJ and ( |92] | are violated, and steady approximation is not valid for radii of the order of Rg-. 
Even in the very beginning of bubble regular growth, when R = IR^, the steady condition ( |96] l is fulfilled only at its 
breaking point: the values of ^'^^ exceed the value of i'^^, but these values are of the same order of magnitude. 
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10. Conclusions 

In the presented paper we obtained the equations for the bubble growth dynamics in the gas solution with Henry's 
and Sievert's solubility laws. We solved these equations analytically for case of bubble growth in strongly super- 
saturated solution. The equation for the Sievert's law was solved also for the case of bubble dissolution in the pure 
solvent. 

We showed that, irrespective of the to gas solubility law, three characteristic stages could be marked out in the 
growth dynamics. During the first stage the bubble radius is growing with an increasing rate. On the second stage the 
growth rate decreases. The third stage, when the growth rate continues to decrease, begins when the Laplace pressure 
inside the bubble becomes comparable with the external pressure of the solution. We demonstrated that during the 
first two stages the time dependence of the bubble radius is different for the cases of Henry's and Sievert's laws, while 
during the third stage this distinction is no longer observed. For both Henry's and Sievert's laws we obtained intervals 
within which the bubble radius changes on each stage, as well as durations of consecutive stages. 

While obtaining the dynamic equations we assumed the diffusion flux to be steady. We obtained conditions when 
this steady approximation is applicable. We showed that usually, as the radius of the bubble increases, the steady 
regime of bubble growth gradually gives way to the nonsteady one. Application of the obtained conditions for the 
volcanic system consisting of water vapor dissolved in the silicate melt showed that the process in such system as a 
rule cannot be considered as steady. 



Appendix A: Effect of solvent viscosity on the bubble growth 

In the current paper we neglected the solvent viscosity. Its influence on bubble dynamics can be estimated using 
Rayleigh-Plesset equation (see e.g. IfTTll ). To take the solvent viscosity into account one needs to replace Eq. (jTjl with 
the following equality: 

2o- R 

P« = n+ — +4/7-, (103) 

where rj is the dynamic viscosity of the solvent. The inertial terms in Rayleigh-Plesset equation are negligible for any 
reasonable bubble growth rate. 

To neglect the viscous term in Eq. ( I103l l (the third addend) in comparison with the surface tension term (the second 
addend) the following strong inequality has to be fulfilled: 

77 « ^ . (104) 

This inequality, evidently, is equivalent to inequaUty f e » 1, where Pe is Peclet number 

The higher the bubble growth rate R is, the stronger the inequality (1104b is. The strongest condition takes place 
when R = 7?,„. Using Eqs. (|49] | and dSTl l in Eqs. (l44l i and (|52] | correspondingly, for both Henry's and Sievert's laws we 
have 

. I Dsi: 

R<RL,,,,-^- (105) 
Therefore, using Eq. ([13}, we can rewrite strong inequality ( llQ4l t as 

77 <K -, (106) 

where multiplier 2/3 in the r h. s. is omitted for shortness. 

Let us estimate the value in the r. h. s. of inequality (1106b . Typical values of surface tension both for water 
and for volcanic systems ||6l are cr ~ 10"' Nm '; both for Henry's and Sievert's laws ~ 10; diffusion coefficient 
D ~ 10"" m^s"' lO; pressure 11 ~ 10^ Pa. Substituting these values, we have: 

?7<KlO^Pas. (107) 

It should be noted that "common" liquids at normal conditions always satisfy this condition: for water we have 
T] ~ 10^^ Pa s and even for glycerol ?; ~ 1 Pa s ifTSll . For volcanic systems the values of viscosity that satisfy strong 
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inequality ( I107l i are quite typical when Si02 content is not too high (basalt, andesite and dactite melts) 111. However, 
for rhyolite melts (~ 70% SiOa) viscosity can reach the values of 10^ Pa s JS); and, therefore, effect of solvent 
viscosity has to be taken into account. 

We do not discuss here the oscillating settlement of the mechanical equilibrium between the bubble and the 
solution: as it was shown in |fT9l this settlement occurs much faster than the settlement of chemical equilibrium, 
unless the liquid viscosity is extremely low. 



Appendix B: Dissolution of the gas bubble in a pure solvent: Sievert's solubility law 

Eq. ( l23T l allows us to obtain the radius-time relation for the bubble of arbitrary initial size R\,=o = Ri which is put 
in the pure solvent and also the time of its dissolution tj, R\,=,. = 0. Separating variables, we can rewrite Eq. (l2Jt as 



1 + ^ 



RdR = -Dsdt , 



(108) 



or equivalently as 



1 + 1?JL _ i ^ 
2 R 2 



RdR = -Dsdt 



Using the variable x = -^1 + |^ for integrating in Eq. (|109t instead of R we have 



(109) 



zRt 



(x2 - 1)3 (x^ - 1)2 



dx - Dsdt . 



(110) 



Eq. ( II 101 ) can be easily integrated and therefore we obtain 

1 I 3 1 I 3 3 

-R Jr^ + -R^R - -R^ Jr^ + -R^R rI In 

2V 2 8 V 2 32 



4- R 4 1 / , 3 

Jr^ + -R^R + 1 

3R,r 3R^y 2 ^ 



= Dst 



(111) 



Eq. (Ill II) gives us an explicit relation between the bubble radius R and time t. 

Let us find the time t^ of total dissolution {R — > 0) from arbitrary initial radius /?, 



td 



1 

IJs 



1 



1 



Ri JR^ + -R^Ri - -Rir \ R + -RcrRi ^^In 

2 y ' 2 8 \' 2 32 



4 Ri 4 1 / , 3 

Jr^ + -Ro-Ri + 1 

3R^ 3R^y ' 2 " ' 



(112) 



The latter expression can be simplified for the two particular cases and short analytical expressions for time t^ can be 
obtained. The first case is /?, » Ry. 

Rj 

td^WTT- (113) 



And the second is Ri <s; Ra-: 



2Ds 



3/2 p3/2pl/2 



Ds 



Let us also obtain the expression for dissolution time from = Rg- 

td = 0.435 . 
Ds 



(114) 



(115) 
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